Differential usage of DNA modifications in neurons, astrocytes, and microglia.

BACKGROUND: Cellular identity is determined partly by cell type-specific epigenomic profiles that regulate gene expression. In neuroscience, there is a pressing need to isolate and characterize the epigenomes of specific CNS cell types in health and disease. In this study, we developed an in vivo tagging mouse model (Camk2a-NuTRAP) for paired isolation of neuronal DNA and RNA without cell sorting and then used this model to assess epigenomic regulation, DNA modifications in particular, of gene expression between neurons and glia. RESULTS: After validating the cell-specificity of the Camk2a-NuTRAP model, we performed TRAP-RNA-Seq and INTACT-whole genome oxidative bisulfite sequencing (WGoxBS) to assess the neuronal translatome and epigenome in the hippocampus of young mice (4 months old). WGoxBS findings were validated with enzymatic methyl-Seq (EM-Seq) and nanopore sequencing. Comparing neuronal data to microglial and astrocytic data from NuTRAP models, microglia had the highest global mCG levels followed by astrocytes and then neurons, with the opposite pattern observed for hmCG and mCH. Differentially modified regions between cell types were predominantly found within gene bodies and distal intergenic regions, rather than proximal promoters. Across cell types there was a negative correlation between DNA modifications (mCG, mCH, hmCG) and gene expression at proximal promoters. In contrast, a negative correlation of gene body mCG and a positive relationship between distal promoter and gene body hmCG with gene expression was observed. Furthermore, we identified a neuron-specific inverse relationship between mCH and gene expression across promoter and gene body regions. CONCLUSIONS: Neurons, astrocytes, and microglia demonstrate different genome-wide levels of mCG, hmCG, and mCH that are reproducible across analytical methods. However, modification-gene expression relationships are conserved across cell types. Enrichment of differential modifications across cell types in gene bodies and distal regulatory elements, but not proximal promoters, highlights epigenomic patterning in these regions as potentially greater determinants of cell identity. These findings also demonstrate the importance of differentiating between mC and hmC in neuroepigenomic analyses, as up to 30% of what is conventionally interpreted as mCG can be hmCG, which often has a different relationship to gene expression than mCG.

Differential usage of DNA modifications in neurons, astrocytes, and microglia.

Background: Cellular identity is determined partly by cell type-specific epigenomic profiles that regulate gene expression. In neuroscience, there is a pressing need to isolate and characterize the epigenomes of specific CNS cell types in health and disease. This is especially true as for DNA modifications where most data are derived from bisulfite sequencing that cannot differentiate between DNA methylation and hydroxymethylation. In this study, we developed an in vivo tagging mouse model (Camk2a-NuTRAP) for paired isolation of neuronal DNA and RNA without cell sorting and then used this model to assess epigenomic regulation of gene expression between neurons and glia. Results: After validating the cell-specificity of the Camk2a-NuTRAP model, we performed TRAP-RNA-Seq and INTACT whole genome oxidative bisulfite sequencing to assess the neuronal translatome and epigenome in the hippocampus of young mice (3 months old). These data were then compared to microglial and astrocytic data from NuTRAP models. When comparing the different cell types, microglia had the highest global mCG levels followed by astrocytes and then neurons, with the opposite pattern observed for hmCG and mCH. Differentially modified regions between cell types were predominantly found within gene bodies and distal intergenic regions, with limited differences occurring within proximal promoters. Across cell types there was a negative correlation between DNA modifications (mCG, mCH, hmCG) and gene expression at proximal promoters. In contrast, a negative correlation of mCG with gene expression within the gene body while a positive relationship between distal promoter and gene body hmCG and gene expression was observed. Furthermore, we identified a neuron-specific inverse relationship between mCH and gene expression across promoter and gene body regions. Conclusions: In this study, we identified differential usage of DNA modifications across CNS cell types, and assessed the relationship between DNA modifications and gene expression in neurons and glia. Despite having different global levels, the general modification-gene expression relationship was conserved across cell types. The enrichment of differential modifications in gene bodies and distal regulatory elements, but not proximal promoters, across cell types highlights epigenomic patterning in these regions as potentially greater determinants of cell identity.

Differential Regulation of Mouse Hippocampal Gene Expression Sex Differences by Chromosomal Content and Gonadal Sex.

Common neurological disorders, like Alzheimer's disease (AD), multiple sclerosis (MS), and autism, display profound sex differences in prevalence and clinical presentation. However, sex differences in the brain with health and disease are often overlooked in experimental models. Sex effects originate, directly or indirectly, from hormonal or sex chromosomal mechanisms. To delineate the contributions of genetic sex (XX v. XY) versus gonadal sex (ovaries v. testes) to the epigenomic regulation of hippocampal sex differences, we used the Four Core Genotypes (FCG) mouse model which uncouples chromosomal and gonadal sex. Transcriptomic and epigenomic analyses of ~ 12-month-old FCG mouse hippocampus, revealed genomic context-specific regulatory effects of genotypic and gonadal sex on X- and autosome-encoded gene expression and DNA modification patterns. X-chromosomal epigenomic patterns, classically associated with X-inactivation, were established almost entirely by genotypic sex, independent of gonadal sex. Differences in X-chromosome methylation were primarily localized to gene regulatory regions including promoters, CpG islands, CTCF binding sites, and active/poised chromatin, with an inverse relationship between methylation and gene expression. Autosomal gene expression demonstrated regulation by both genotypic and gonadal sex, particularly in immune processes. These data demonstrate an important regulatory role of sex chromosomes, independent of gonadal sex, on sex-biased hippocampal transcriptomic and epigenomic profiles. Future studies will need to further interrogate specific CNS cell types, identify the mechanisms by which sex chromosomes regulate autosomes, and differentiate organizational from activational hormonal effects.

Inducible cell-specific mouse models for paired epigenetic and transcriptomic studies of microglia and astroglia.

Epigenetic regulation of gene expression occurs in a cell type-specific manner. Current cell-type specific neuroepigenetic studies rely on cell sorting methods that can alter cell phenotype and introduce potential confounds. Here we demonstrate and validate a Nuclear Tagging and Translating Ribosome Affinity Purification (NuTRAP) approach for temporally controlled labeling and isolation of ribosomes and nuclei, and thus RNA and DNA, from specific central nervous system cell types. Analysis of gene expression and DNA modifications in astrocytes or microglia from the same animal demonstrates differential usage of DNA methylation and hydroxymethylation in CpG and non-CpG contexts that corresponds to cell type-specific gene expression. Application of this approach in LPS treated mice uncovers microglia-specific transcriptome and epigenome changes in inflammatory pathways that cannot be detected with tissue-level analysis. The NuTRAP model and the validation approaches presented can be applied to any brain cell type for which a cell type-specific cre is available.

Early-life DNA methylation profiles are indicative of age-related transcriptome changes.

BACKGROUND: Alterations to cellular and molecular programs with brain aging result in cognitive impairment and susceptibility to neurodegenerative disease. Changes in DNA methylation patterns, an epigenetic modification required for various CNS functions are observed with brain aging and can be prevented by anti-aging interventions, but the relationship of altered methylation to gene expression is poorly understood. RESULTS: Paired analysis of the hippocampal methylome and transcriptome with aging of male and female mice demonstrates that age-related differences in methylation and gene expression are anti-correlated within gene bodies and enhancers. Altered promoter methylation with aging was found to be generally un-related to altered gene expression. A more striking relationship was found between methylation levels at young age and differential gene expression with aging. Highly methylated gene bodies and promoters in early life were associated with age-related increases in gene expression even in the absence of significant methylation changes with aging. As well, low levels of methylation in early life were correlated to decreased expression with aging. This relationship was also observed in genes altered in two mouse Alzheimer's models. CONCLUSION: DNA methylation patterns established in youth, in combination with other epigenetic marks, were able to accurately predict changes in transcript trajectories with aging. These findings are consistent with the developmental origins of disease hypothesis and indicate that epigenetic variability in early life may explain differences in aging trajectories and age-related disease.

Tamoxifen induction of Cre recombinase does not cause long-lasting or sexually divergent responses in the CNS epigenome or transcriptome: implications for the design of aging studies.

The systemic delivery of tamoxifen (Tam) to activate inducible CreERT2-loxP transgenic mouse systems is now widely used in neuroscience studies. This critical technological advancement allows temporal control of DNA-cre recombination, avoidance of embryonically lethal phenotypes, and minimization of residual cell labeling encountered in constitutively active drivers. Despite its advantages, the use of Tam has the potential to cause long-lasting, uncharacterized side effects on the transcriptome and epigenome in the CNS, given its mixed estrogen receptor (ER) agonist/antagonist actions. With the welcome focus on including both sexes in biomedical studies and efforts to understand sex differences, Tam administration could also cause sexually divergent responses that would confound studies. To examine these issues, epigenetic and transcriptomic profiles were compared in C57BL/6 J female and male hippocampus, cortex, and retina 1 month after a 5-day Tam treatment typical for cre induction, or vehicle control (sunflower seed oil). Cytosine methylation and hydroxymethylation levels, in both CG and non-CG contexts, were unchanged as determined by oxidative bisulfite sequencing. Long-lasting Tam transcriptomic effects were also not evident/minimal. Furthermore, there is no evidence of sexually divergent responses with Tam administration and Tam did not alter sex differences evident in controls. Combined with recently reported data that Tam alone does not cause long-lasting changes in behavior and neurogenesis, our findings provide confidence that Tam can be used as a cre-recombinase inducer without introducing significant confounds in transcriptomic and epigenomic neuroscience studies, particularly those focused on genomic and transcriptomic aspects of the aging brain.

Exposure to environmental enrichment attenuates addiction-like behavior and alters molecular effects of heroin self-administration in rats.

Environmental factors profoundly affect the addictive potential of drugs of abuse and may also modulate the neuro-anatomical/neuro-chemical impacts of uncontrolled drug use and relapse propensity. This study examined the impact of environmental enrichment on heroin self-administration, addiction-related behaviors, and molecular processes proposed to underlie these behaviors. Male Sprague-Dawley rats in standard and enriched housing conditions intravenously self-administered similar amounts of heroin over 14 days. However, environmental enrichment attenuated progressive ratio, extinction, and reinstatement session responding after 14 days of enforced abstinence. Molecular mechanisms, namely DNA methylation and gene expression, are proposed to underlie abstinence-persistent behaviors. A global reduction in methylation is reported to coincide with addiction, but no differences in total genomic methylation or repeat element methylation were observed in CpG or non-CpG (CH) contexts across the mesolimbic circuitry as assessed by multiple methods including whole genome bisulfite sequencing. Immediate early gene expression associated with drug seeking, taking, and abstinence also were examined. EGR1 and EGR2 were suppressed in mesolimbic regions with heroin-taking and environmental enrichment. Site-specific methylation analysis of EGR1 and EGR2 promoter regions using bisulfite amplicon sequencing (BSAS) revealed hypo-methylation in the EGR2 promoter region and EGR1 intragenic CpG sites with heroin-taking and environmental enrichment that was associated with decreased mRNA expression. Taken together, these findings illuminate the impact of drug taking and environment on the epigenome in a locus and gene-specific manner and highlight the need for positive, alternative rewards in the treatment and prevention of drug addiction.

Caloric restriction mitigates age-associated hippocampal differential CG and non-CG methylation.

Brain aging is marked by cognitive decline and susceptibility to neurodegeneration. Calorie restriction (CR) increases neurogenesis, improves memory function, and protects from age-associated neurological disorders. Epigenetic mechanisms, including DNA methylation, are vital to normal central nervous system cellular and memory functions and are dysregulated with aging. The beneficial effects of CR have been proposed to work through epigenetic processes, but this is largely unexplored. We therefore tested whether life long CR prevents age-related hippocampal DNA methylation changes. Hippocampal DNA from young (3 months) and old (24 months) male mice fed ad libitum and 24-month-old mice fed a 40% calorie-restricted diet from 3 months of age were examined by genome-wide bisulfite sequencing to measure methylation with base specificity. Over 27 million CG and CH (non-CG) sites were examined. Of the ∼40,000 differentially methylated CG and ∼80,000 CH sites with aging, >1/3 were prevented by CR and were found across genomic regulatory regions and gene pathways. CR also caused alterations to CG and CH methylation at sites not differentially methylated with aging, and these CR-specific changes demonstrated a different pattern of regulatory element and gene pathway enrichment than those affected by aging. CR-specific DNA methyltransferase 1 and Tet methylcytosine dioxygenase 3 promoter hypermethylation corresponded to reduced gene expression. These findings demonstrate that CR attenuates age-related CG and CH hippocampal methylation changes, in combination with CR-specific methylation that may also contribute to the neuroprotective effects of CR. The prevention of age-related methylation alterations is also consistent with the prolongevity effects of CR working through an epigenetic mechanism.

Analysis of DNA modifications in aging research.

As geroscience research extends into the role of epigenetics in aging and age-related disease, researchers are being confronted with unfamiliar molecular techniques and data analysis methods that can be difficult to integrate into their work. In this review, we focus on the analysis of DNA modifications, namely cytosine methylation and hydroxymethylation, through next-generation sequencing methods. While older techniques for modification analysis performed relative quantitation across regions of the genome or examined average genome levels, these analyses lack the desired specificity, rigor, and genomic coverage to firmly establish the nature of genomic methylation patterns and their response to aging. With recent methodological advances, such as whole genome bisulfite sequencing (WGBS), bisulfite oligonucleotide capture sequencing (BOCS), and bisulfite amplicon sequencing (BSAS), cytosine modifications can now be readily analyzed with base-specific, absolute quantitation at both cytosine-guanine dinucleotide (CG) and non-CG sites throughout the genome or within specific regions of interest by next-generation sequencing. Additional advances, such as oxidative bisulfite conversion to differentiate methylation from hydroxymethylation and analysis of limited input/single-cells, have great promise for continuing to expand epigenomic capabilities. This review provides a background on DNA modifications, the current state-of-the-art for sequencing methods, bioinformatics tools for converting these large data sets into biological insights, and perspectives on future directions for the field.

Sexually divergent DNA methylation patterns with hippocampal aging.

DNA methylation is a central regulator of genome function, and altered methylation patterns are indicative of biological aging and mortality. Age-related cellular, biochemical, and molecular changes in the hippocampus lead to cognitive impairments and greater vulnerability to neurodegenerative disease that varies between the sexes. The role of hippocampal epigenomic changes with aging in these processes is unknown as no genome-wide analyses of age-related methylation changes have considered the factor of sex in a controlled animal model. High-depth, genome-wide bisulfite sequencing of young (3 month) and old (24 month) male and female mouse hippocampus revealed that while total genomic methylation amounts did not change with aging, specific sites in CG and non-CG (CH) contexts demonstrated age-related increases or decreases in methylation that were predominantly sexually divergent. Differential methylation with age for both CG and CH sites was enriched in intergenic and intronic regions and under-represented in promoters, CG islands, and specific enhancer regions in both sexes, suggesting that certain genomic elements are especially labile with aging, even if the exact genomic loci altered are predominantly sex-specific. Lifelong sex differences in autosomal methylation at CG and CH sites were also observed. The lack of genome-wide hypomethylation, sexually divergent aging response, and autosomal sex differences at CG sites was confirmed in human data. These data reveal sex as a previously unappreciated central factor of hippocampal epigenomic changes with aging. In total, these data demonstrate an intricate regulation of DNA methylation with aging by sex, cytosine context, genomic location, and methylation level.

CNS-wide Sexually Dimorphic Induction of the Major Histocompatibility Complex 1 Pathway With Aging.

The major histocompatibility complex I (MHCI) pathway, which canonically functions in innate immune viral antigen presentation and detection, is functionally pleiotropic in the central nervous system (CNS). Alternative roles include developmental synapse pruning, regulation of synaptic plasticity, and inhibition of neuronal insulin signaling; all processes altered during brain aging. Upregulation of MHCI components with aging has been reported; however, no systematic examination of MHCI cellular localization, expression, and regulation across CNS regions, life span, and sexes has been reported. In the mouse, MHCI is expressed by neurons and microglia, and MHCI components and receptors (H2-K1, H2-D1, ß2M, Lilrb3, Klra2, CD247) display markedly different expression profiles across the hippocampus, cortex, cerebellum, brainstem, and retina. MHCI components, receptors, associated inflammatory transcripts (IL1α, IL1ß, IL6, TNFα), and TAP (transporter associated with antigen processing) components are induced with aging and to a greater degree in female than male mice across CNS regions. H2-K1 and H2-D1 expression is associated with differential CG and non-CG promoter methylation across CNS regions, ages, and between sexes, and concomitant increased expression of proinflammatory genes. Meta-analysis of human brain aging data also demonstrates age-related increases in MHCI. Induction of MHCI signaling could contribute to altered synapse regulation and impaired synaptic plasticity with aging.

Absence of genomic hypomethylation or regulation of cytosine-modifying enzymes with aging in male and female mice.

BACKGROUND: Changes to the epigenome with aging, and DNA modifications in particular, have been proposed as a central regulator of the aging process, a predictor of mortality, and a contributor to the pathogenesis of age-related diseases. In the central nervous system, control of learning and memory, neurogenesis, and plasticity require changes in cytosine methylation and hydroxymethylation. Although genome-wide decreases in methylation with aging are often reported as scientific dogma, primary research reports describe decreases, increases, or lack of change in methylation and hydroxymethylation and their principle regulators, DNA methyltransferases and ten-eleven translocation dioxygenases in the hippocampus. Furthermore, existing data are limited to only male animals. RESULTS: Through examination of the hippocampus in young, adult, and old male and female mice by antibody-based, pyrosequencing, and whole-genome oxidative bisulfite sequencing methods, we provide compelling evidence that contradicts the genomic hypomethylation theory of aging. We also demonstrate that expression of DNA methyltransferases and ten-eleven translocation dioxygenases is not differentially regulated with aging or between the sexes, including the proposed cognitive aging regulator DNMT3a2. Using oxidative bisulfite sequencing that discriminates methylation from hydroxymethylation and by cytosine (CG and non-CG) context, we observe sex differences in average CG methylation and hydroxymethylation of the X chromosome, and small age-related differences in hydroxymethylation of CG island shores and shelves, and methylation of promoter regions. CONCLUSION: These findings clarify a long-standing misconception of the epigenomic response to aging and demonstrate the need for studies of base-specific methylation and hydroxymethylation with aging in both sexes.

Bisulfite oligonucleotide-capture sequencing for targeted base- and strand-specific absolute 5-methylcytosine quantitation.

Epigenetic regulation through DNA methylation (5mC) plays an important role in development, aging, and a variety of diseases. Genome-wide studies of base- and strand-specific 5mC are limited by the extensive sequencing required. Targeting bisulfite sequencing to specific genomic regions through sequence capture with complimentary oligonucleotide probes retains the advantages of bisulfite sequencing while focusing sequencing reads on regions of interest, enables analysis of more samples by decreasing the amount of sequence required per sample, and provides base- and strand-specific absolute quantitation of CG and non-CG methylation levels. As an example, an oligonucleotide capture set to interrogate 5mC levels in all rat RefSeq gene promoter regions (18,814) and CG islands, shores, and shelves (18,411) was generated. Validation using whole-genome methylation standards and biological samples demonstrates enrichment of the targeted regions and accurate base-specific quantitation of CG and non-CG methylation for both forward and reverse genomic strands. A total of 170 Mb of the rat genome is covered including 6.6 million CGs and over 67 million non-CG sites, while reducing the amount of sequencing required by ~85 % as compared to existing whole-genome sequencing methods. This oligonucleotide capture targeting approach and quantitative validation workflow can also be applied to any genome of interest.

Targeted DNA methylation analysis by next-generation sequencing.

The role of epigenetic processes in the control of gene expression has been known for a number of years. DNA methylation at cytosine residues is of particular interest for epigenetic studies as it has been demonstrated to be both a long lasting and a dynamic regulator of gene expression. Efforts to examine epigenetic changes in health and disease have been hindered by the lack of high-throughput, quantitatively accurate methods. With the advent and popularization of next-generation sequencing (NGS) technologies, these tools are now being applied to epigenomics in addition to existing genomic and transcriptomic methodologies. For epigenetic investigations of cytosine methylation where regions of interest, such as specific gene promoters or CpG islands, have been identified and there is a need to examine significant numbers of samples with high quantitative accuracy, we have developed a method called Bisulfite Amplicon Sequencing (BSAS). This method combines bisulfite conversion with targeted amplification of regions of interest, transposome-mediated library construction and benchtop NGS. BSAS offers a rapid and efficient method for analysis of up to 10 kb of targeted regions in up to 96 samples at a time that can be performed by most research groups with basic molecular biology skills. The results provide absolute quantitation of cytosine methylation with base specificity. BSAS can be applied to any genomic region from any DNA source. This method is useful for hypothesis testing studies of target regions of interest as well as confirmation of regions identified in genome-wide methylation analyses such as whole genome bisulfite sequencing, reduced representation bisulfite sequencing, and methylated DNA immunoprecipitation sequencing.